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Abstract 

Pairing transition at finite temperature was investigated by the shell model and 
BCS calculations. The definitive signature of pairing transition is identified by a 
"transition temperature" Tf estimated from a "thermal" odd-even mass difference, 
while there is no sharp phase transition because of the finiteness of nucleus. It is 
found that Tt is in good agreement with predictions of critical temperature Tc in 
the BCS approximation, and the pairing correlations almost vanish at two points 
of the transition temperature T ss 2Tt. The BCS calculations show that the critical 
temperature Tc increases with increasing deformation. 
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1 Introduction 

Pairing correlations are one of the fundamental properties of nuclei. The odd- 
even mass difference observed in nuclear masses, and the energy gap of even- 
even nuclei are well known as signatures of pairing effects [1]. The Bardeen- 
Cooper-Schriffer (BCS) theory [2] of superconductivity has been applied to 
such nuclear problems at zero temperature [3]. The thermodynamical prop- 
erties of nuclear pairing were investigated using the BCS theory in the study 
of hot nuclei [4,5]. Breaking of the Cooper pairs is expected to occur above a 
certain critical temperature. Infinite Fermi systems show a sharp phase tran- 
sition from a super-fluid phase to a normal-fluid one at a critical temperature 
Tc, where the pairing gap vanishes and the heat capacity exhibits a singularity. 
For metal superconductors, the critical temperature is Tc ~ 0.57A MeV [2], 
where the pairing gap A is calculated using the BCS equation. In this case, we 
can image a well deformed nucleus for which level density is almost uniform 
and average level spacing is very small compared with the pairing gap A 
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1 MeV for the rare earth nuclei) [6]. On the other hand, the pairing correla- 
tions vanish at Tc ~ 0.50A MeV [7] in a simple pairing model with half-filled 
degenerate shells. 

However in the case of the finite Fermi system like a nucleus, since the nuclear 
radius is much smaller than the coherence length of the Cooper pair, statistical 
fluctuations beyond the mean field become larger. The fluctuations smooth out 
the sharp phase transition, and then the pairing correlations do not vanish 
but decrease with increasing temperature. There are many approaches for 
treating the fluctuations beyond the mean field. A signature of the pairing 
transition at finite temperature might be a peak of the heat capacity as a 
function of temperature [8]. In fact, it has recently been reported [9,10] that 
the canonical heat capacities extracted from observed level densities in ^^^Dy, 
^^^Er and ^'''^Yb form an S shape around T 0.5 MeV, which is interpreted 
as the breaking of nucleon Cooper pairs and the pairing transition because 
the critical temperature corresponds to Tc ~ 0.57A ^ 0.5 MeV in the BCS 
theory. The odd-even effect that shows that the heat capacity of an odd-mass 
nucleus is smaller than those of the adjacent even- mass nuclei can be found 
around the critical temperature Tc- 

The spherical shell model approaches could be more appropriate for describ- 
ing various aspects of nuclear structure. The large fluctuations can be taken 
into account beyond the mean fleld in the shell model calculations. For the de- 
scription of nuclear properties at higher temperatures, one needs a large model 
space. Recently the shell model Monte Carlo (SMMC) calculations have been 
performed with the pure pairing force [11] in a large model space, and with 
the pairing and multipole-multipolc forces [12] in the fp + gg/2 shell model 
space for the even- and odd-mass Fc isotopes. It was shown that the pairing 
correlations would be important only at low temperatures and at low exci- 
tation energies [11], and that the suppression of pairing correlations due to 
finite temperature appears as the S shape of the heat capacity around the 
temperature T ~ 0.8 MeV [12]. The model space was quite recently extended 
to examine the partition functions and the level densities up to the higher 
temperature T = 4.0 MeV using an independent-particle approximation [13] 
combined with the SMMC method. 

This paper is organized as follows. In Section 2, we carry out the shell model 
calculations in the large model space {sd + fp + 51/2^^5/2) using a spherical 
Woods-Saxon potential [14], where we adopt the independent-particle approx- 
imation [13] and combine it with the shell model calculations in the sd-shell, 
though there are no contributions from many-body correlations in the fp shell 
or from the coupling between the fp and sd shells. The pairing correlations 
within the heat capacity or "thermal" odd-even mass difference as a func- 
tion of flnite temperature are estimated, and a "transition temperature" Tt 
is introduced to explain the pairing transition. Then the definitive signatures 
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of the pairing transition are identified. The critical temperature Tc is also 
evaluated using the BCS calculations with the axially deformed Woods-Saxon 
potential and a pairing residual interaction in Section 3. We discuss the re- 
lations between the transition temperature Tt and the critical temperature 
Tc- The systematic behavior of the critical temperature Tc is examined over 
a wide range of nuclei. We discuss the effects of deformation on the critical 
temperature. Concluding remarks are given in Section 4. 



2 Pairing transition in shell model calculations 



We start from the canonical partition function defined by 



Z{T) = Tr(e-^/^) = e"^^/^, (1) 

i=0 

where Ei is the energy of the ith eigenstate for the Hamiltonian if of a system. 

The large matrix of the Hamiltonian H is diagonalizcd to obtain all the eigen- 
values -Ej, and the partition function in the canonical ensemble is calculated 
from Eq. (1). Then, any thermodynamical quantities 0{T) can be evaluated 
from 

0(T) = (O) = Tr(Oe-^/^)/Z(T), (2) 

where (O) stands for the average value of operator O over the range of eigen- 
states. For instance, the thermal energy is expressed as 



E{Z,N,T)^{H)^Y. E.e-'^^'^/ZiT) . (3) 
The heat capacity is then given by 



We now introduce the thermal odd-even mass difference for neutrons defined 
by the following three-point indicator: 



A(f) (Z, TV, T) = ^-^[E{Z, TV + 1, T) 



-2E{Z,N,T) + E{Z,N -l.T)]. (5) 
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Fig. 1. Heat capacity (a) and thermal odd-even mass difference (b) as a function of 
temperature. The sohd circles denote the values for even-even nucleus ^^Mg. The 
open circles denote the average values of the neighboring odd nuclei ^^Mg and ^^Mg. 

The odd-even mass difference at zero temperature is known theoretically and 
experimentally as an important quantity in evaluation of the pairing correla- 
tions in a nucleus. The thermal odd-even mass difference is also an indicator 
of the pairing correlations at finite temperature, and can be obtained from the 
experimental energies and the level density as well as the heat capacity. 



We now calculate the heat capacity (4) and the thermal odd-even mass differ- 
ence (5) for ^°^^^Ne and ^^^^^Mg in the large model space {sd + fp + 51/2(^5/2) 
using an independent-particle approximation [13]. The Woods-Saxon potential 
plus a spin-orbit interaction is diagonalized in a basis of harmonic-oscillator 
(H.O.) eigenf unctions, and then the single-particle energies are obtained [14]. 
The Woods-Saxon parameters are chosen so as to reproduce the single-particle 
energies estimated from ^^O. A number of unbound states {E > 0) are ob- 
tained due to the expansion in a finite number of H.O. eigenfunctions. The 
resonances with narrow widths are important for neutron-rich nuclei such as 
^^Cr. It was shown [13], however, that while the narrow resonances make an 
important contribution to the partition function, their width can be ignored. 
In this paper, we ignore the widths of resonances and the continuum states. 
We carried out the shell model calculations in the sii-shell model space using 
the USD interaction [15], and calculated the correlated thermal energy E^^t^ 
using Eq. (3). An enhancement of the heat capacity was found in even Mg 
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nuclei around T ~ 1.5 MeV. This enhancement is interpreted as a reduction 
in the pairing correlations. However, because the calculation was restricted to 
a finite space {sd shell), the calculated heat capacity reached its maximum 
around temperatures of ~ 1.5 MeV. Therefore, we extended the model space 
to sd+fp+si/2d5/2 to display a much broader temperature range (~ 3.5 MeV). 
We combine the small model space sd of interaction effects with a large-space 
calculation of the independent-particle thermal energy Esp as follows [13]: 

E — E^,tr + Egp — Egp^trj (6) 

where Esp^tr is the independent-particle thermal energy in the sd space. Then 
the heat capacity and the thermal odd-even mass difference are obtained 
using Eqs. (4) and (5), respectively. As typical examples, the heat capac- 
ity and the thermal odd-even mass difference in ^'''"^^Mg are plotted in Fig. 
1 (a) and (b), where the average values of ^'''Mg and ^^Mg are shown for 
odd-mass nuclei. The heat capacity C(^^Mg) is larger than the average value 
C(27'29Mg) = [C(^^Mg) + C(29Mg)]/2 in T = 0.8 - 3.0 MeV, although it does 
not reach a peak. In other nuclei, we also found an odd-even effect where 
the average heat capacity of odd-mass nuclei is significantly lower than that 
of the adjacent even-mass nucleus. In Fig. 1 (a), one can see that with the 
usual relation C = 2aT the parameter a has roughly the empirical value of 
A/8 MeV^^ [16] for light nuclei, and is considerably larger than the Fermi- 
gas model value A/15 MeV^^. Figure 1 (b) shows the thermal odd-even mass 
difference defined by Eq. (5) for ^^Mg and ^'''■^^Mg. We find a gradual de- 
crease of the thermal odd-even mass difference as a function of temperature, 
which is interpreted as a gradual breaking of nucleon Cooper pairs and the 
decline of pairing correlations. Figure 1 (b) also displays that A^f)(^^Mg) is 
larger than A(f)(2'^'29Mg) = [A'^^'^'^Mg) + A(f)(29Mg)]/2. At zero temper- 
ature, this is well known as the odd-even staggering of binding energies, 
which refiects the stronger binding of even-particle-number systems than the 
odd-particle-number neighbors [17,18]. Their analyses demonstrated that the 
odd-even mass difference for odd-particle number is an excellent measure of 
pairing correlations, although it is still controversial [19]. The symmetric fil- 
ter Se = 2A(f)(2m) - A(f)(2m - 1) - A(f)(2m + 1) extracts the effective 
single-particle spacing from the measured binding energies of deformed nuclei. 
We can see that the difference, 6e/2, between A[f\'^^Mg) and A^f\'^'^'^^Mg) 
decreases gradually as temperature increases. This means that the effective 
single-particle spacing decreases as temperature increases. As one can see in 
Fig. 1 (b), A(f)(27.29Mg) is almost zero aX T ^ 3.0 MeV, and the pairing 
correlations vanish at two points on the transition temperature T^. On the 
other hand, A^/^\^^Mg) still remains at this temperature. We now identify 
an inffection point of the curve A^^-'('^''''^'^Mg) in Fig. 1 (b) gnature of 

pairing transition, and call it "transition temperature" T^. To see more precise 
position of the inflection point, we differentiate A^^^ (^'''■^^Mg) with respect to 
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Fig. 2. Derivative of the thermal odd-even mass difference. 

temperature T. As seen in Fig. 2, — 9A^^) (^^'^^Mg)/(?T has a peak at tem- 
perature Tc ~ 1.3 MeV corresponding to the transition temperature T^, and 
the transition temperature Tt for ^^'^^Mg is quite close to that for ^^Mg. It is 
very important to note that the difference between the two curves of the heat 
capacities in Fig. 1 (a) is equal to — SA^f ^ (^^Mg) / dT 



f)Ai^)(y AT T] 1 
dT = i-^f{CiZ, N, T) - -[C{Z, N + 1,T) 

+CiZ,N-l,T)]}. (7) 

Thus, the thermal odd-cvcn mass difference is a good indicator for the pairing 
transition at the transition temperature estimated from a peak of — SA^f^ /dT. 
It is important to note that the position of the peak does not change upon 
extension of the model space in the independent-particle approximation, which 
is different from the result of Ref. [11]. 



3 BCS calculations 



In infinite systems, the BCS critical temperature is known to be proportional 
to the pairing gap Tc — 0.57A. However, the relation is not axiomatic in a finite 
system like a nucleus. It seems that the peaks of the S-shaped heat capacities 
in ^^^Dy, ^^^Er, and ^^^Yb correspond to the critical temperature = 0.57A 
in the BCS prediction for infinite systems. It is therefore interesting to examine 
the correspondence using the BCS calculations over a wide range of nuclei. In 
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the BCS approximation, the pair gap A at finite temperature is obtained by 
using the following gap equation 



GEV^, (8) 



where are the single-particle energies, = y {et — + the quasiparti- 
cle energies, and fk = {l + e^'^^'^)~^ the Fermi-Dirac quasiparticle occupancies. 
The chemical potential /i is determined by the number constraint 



^-E[l- V^(l-2/.)], (9) 
fe>0 ^'^ 

where the pairing force strength G is chosen so as to reproduce the experi- 
mental odd-even mass difference at zero temperature. In this calculation, we 
use the single-particle energies extracted from an axially deformed Woods- 
Saxon potential with spin-orbit interaction [14]. We chose the parameters 
so as to accommodate experimental singlc-particlc energies extracted from 
energy levels of odd nucleus with double-closed core plus one neutron, i.e., 
^"^Ca, ^°^Sn, and ^^'^Pb. The deformation takes into account the effects of a 
quadrupole-quadrupole interaction in the mean-field approximation. The de- 
formation parameter in even-even nuclei can be estimated from B(E2) — 
[{'i/ATx)Zer'^A^^^(3]^ using experimental B{E2) values. If they are not avail- 
able, we adopt the empirical formula B{E2) ^ 3.27E~+^Z'^A~^-^^ for the 
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Fig. 3. Thermal odd-even mass difference. The solid circles represent the even-even 
nucleus ^^Mg, and the dotted line the pairing gap in the BCS approximation. 
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Fig. 4. Transition temperature Tt (solid symbols) and critical temperature Tc (open 
symbols) as a function of neutron number in the Ne and Mg isotopes. 

B(E2) values where £'2+ is the energy of first excited 2"*" state in even-even 
nucleus. The number of two-fold degenerate active orbitals is chosen as 20 cor- 
responding to N=2,3 shells above ^'^Ca core in sd-nuclei, and 20 corresponding 
to N=3,4 shells above ^''^Sn core in /p-nuclei for each proton and neutron. For 
heavy nuclei, the numbers of active proton and neutron orbitals arc taken as 
30 corresponding to N=4,5 and N=5,6 shells above ^osp]-, ^ore, respectively. 
The BCS pairing gap for ^^Mg is shown in Fig. 3. The value of A^^) decreases 
with increasing temperature, and vanishes at Tc ~ 1.4 MeV for ^^Mg. We can 
see that the critical temperature for ^^Mg is very close to the peak of the heat 
capacity. 

For the Ne and Mg isotopes examined, we examine the temperature about 
the peak of —dA^^/dT in the shell model calculation and the critical tem- 
perature Tc in the BCS. Figure 4 suggests that the critical temperature 
can be identified with the transition temperature Tt even though the pairing 
correlations do not vanish. It is very important to note that the N — Z nuclei 
^°Ne, ^■^Mg show a high transition temperature Tt ^ 2.7 MeV, while the other 
N > Z nuclei have Tt ~ 1.3 MeV. In our previous paper [20], we suggested 
that proton-neutron {p-n) correlations give rise to large odd-even mass differ- 
ence in = Z nuclei. The cooperation of the p-n pairing correlations with the 
like-nucleon correlations would increase the transition temperature as well. 

We have found correspondence between the "transition temperature" and the 
BCS critical temperature in light nuclei. It is worth examining whether or 
not the correspondence obtained in light nuclei also holds in heavy nuclei 
over a wide range. Carrying out exact shell-model calculations is, however. 
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Fig. 5. The upper graph (a) shows the neutron pairing gap A versus temperature 
T. The lower graph (b) shows the neutron pairing gap at zero temperature versus 
the critical temperature Tc in the BCS calculation. 

possible only for light nuclei but not for heavy nuclei. We apply the "cor- 
respondence formula" found in the BCS approximation for the Mg and Ne 

isotopes to heavy nuclei. Since the transition temperature Tf can be regarded 
as the critical temperature T^. from the above discussions, it is now interesting 
to investigate the neutron gap and the critical temperature in the BCS calcu- 
lation solving the gap equation (8) with the number equation (9) over a wide 
range of even-even nuclei, i.e., ^'^'^^O, ^''-^'Ne, ^^~^^Mg, ^""-^^Si, ^^-a^S, ^^-^Ar, 

42-46Ca, ^^-'^^Ti, 48-52Cr, 54-66Fe, 104-124gj^^ lOS-lM^e, 124,130Ce, i32Nd, i34Sm, 

i4o,i44,154qj^ i62j3y^ i66j^^^ ^^^^ i72Yb. The pairing force strength G is deter- 
mined so as to reproduce the experimental odd-even mass difference at zero 
temperature. The neutron pairing gap vanishes at the critical temperature for 
each nuclei as seen in the upper graph (a) of Fig. 5. The critical temperatures 
of ^^^Dy, -"^^^Er, and -"^^^Yb are around Tc ~ 0.5 MeV, which corresponds to the 
peak of the experimental heat capacity [9,10]. In addition, the critical temper- 
ature Tc ~ 0.8 MeV for ®^Fe agrees with the temperature of peak in the heat 
capacity obtained using the SMMC calculations [11,12]. These results support 
the conjecture that the pairing transition temperature T^ corresponds to the 
critical temperature T^. In the lower graph (b) of Fig. 5, the neutron gaps at 
zero temperature are plotted as a function of Tc and the two lines are denoted 
as a guide eye. One is the line ~ 0.57A derived from the BCS theory [2]. 
The other is that of a simple pairing model with half-filled degenerate shells. 
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Fig. 6. Deformation dependence of the critical temperature Tc in the BCS calcula- 
tion. The upper graph (a) shows the neutron gap in /3=0.0, 0.08, and 0.334. The 
lower graph (b) shows the ratio of the critical temperature Tc and the gap A„. 

where the pairing gaps vanish at Tc ~ 0.50 A MeV [7]. Almost all the plots 
lie between these two lines, and those for deformed nuclei are very close to 
the line Tc = 0.57A. For instance, the critical temperature of ^^^Er exhibits 
Tc = 0.58 MeV where A=1.02 MeV and the deformation /3=0.334. Since a well 
deformed nucleus can be regarded as a system with almost uniform level den- 
sity and a small average level spacing compared with the gap A 1 MeV for 
the rare earth nuclei) [6], this result is reasonable. This is simply the case for 
the metal superconductors. The critical temperature is given by Tc ~ 0.57A 
(A = ue^^^^^) where u is the phonon energy and p the average level density at 
the Fermi surface, and Tc is sensitive to the level density. On the other hand, 
the simple pairing model with half-filled degenerate shells leads to Tc ~ 0.50A 
MeV [7]. Thus it seems that the critical temperature depends on the degree of 
deformation. In fact. Fig. 6 shows that the critical temperature Tc increases 
with increasing deformation in ^^^Er. 



4 Concluding remcirks 



In conclusion, we have studied the pairing transition at finite temperature 
using the shell model calculations. We have demonstrated that the thermal 
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odd-even mass difference is a good indicator for tlie pairing transition at finite 
temperature as well as the usual one at zero temperature. We suggest that 
the pairing correlations can be estimated from the measured level densities of 
nuclei with neutron number A^^+l, A^, and A^— 1, for instance, ^'''°Yb, ^^^Yb, and 
^^^Yb. It was shown that the transition temperature Tt corresponding to the 
inflection point of the curve A^^^ is almost identical to the critical temperature 
Tc in the BCS method. The pairing correlations almost vanish at two points 
on the transition temperature T^. The critical temperature Tc depends on 
the deformation of a nucleus, and increases with increasing deformation. The 
critical temperature Tc of deformed nuclei follows Tc ~ 0.57 A MeV as the 
case of the metal superconductors. The transition temperature in the case of 
N = Z nuclei is comparatively higher than those of the neighboring N > Z 
nuclei. The p-n pair correlations seems to contribute to the increase of the 
transition temperature. In addition, as the p-n pairing is crucial for formation 
of an a-like correlated structure in the N — Z nuclei [21], we can expect the 
breaking of p-n pairs as well as like-nucleon pairs at a certain temperature. 
Further studies in this direction are in progress. 
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